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Abstract: In the framework of the QCD light-cone sum rules (LCSRs) we present the 
analysis of all B , B s —> rf'> and D , D s — > rf 1 ' 1 form factors (/ + , /° and f T ) by including m^,) 
corrections in the leading (up to the twist-four) and next-to-leading order (up to the twist- 
three) in QCD, and two-gluon contributions to the form factors at the leading twist. The 
SU(3)-flavour breaking corrections and the axial anomaly contributions to the distribution 
amplitudes are also consistently taken into account. The complete results for the /° and 
f T form factors of B , B s -A 77 W and D , D s —> 77 ^ relevant for processes like B —> t{^tv t or 
b s -+ 7/Wz+r are given for the first time, as well as the two-gluon contribution to the ten¬ 
sor form factors. The values obtained for the f + form factors are as follows: /A (0) = 


o-i68±8;g“, \i 


,+0.036 


/b 3 t/'(0) — 0-252^o;q 20 an d 

113 


B sT? (0)I — 0.212 0 ; 0 13 , /^/(O) — 0. 032 ; — <J " i '-'^- 0.020 


— 0.429lon4i, |/i„(0)| — 0-495 ^o;o29> /d?/(0) — 0.292+2^4, /d s7? '(0) — 

Also phenomenological predictions for semileptonic B, B s —> 77 ^ and D, D s —> 77 ^ decay 
modes are given. 


1 Corresponding author. 
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1 Introduction 

In the view of the numerous precise new measurements of two-body nonleptonic and 

semileptonic B,B S and D,D S decays to r performed by BaBar and Belle recently [1] 

and the upcoming experimental precision in the next-generation experiments it is timely 

to provide precise predictions for B,B S —> rand D,D S —> r/^ 7 -* form factors for analysis 
of these decays. The form factors parametrize hadronic matrix elements of quark currents 

and describe the long-distance QCD effects in semi-leptonic and non-leptonic decays. 

All those decays are important for testing and understanding the Standard Model 
flavour interactions, in particular for our understanding of the QCD dynamics in the flavour 
physics as well as the flavour mixing given by the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix. The B, B s and D, D s decays to r/, rj pseudoscalar mesons can be used to 
shed some light on both of these phenomena. 
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Specially, the decays B,B S —> X c5 P, where X c5 = J/tf>,ip,rf c ,XdO,ci,c 2 ,h c and P is the 
light pseudoscalar meson P = n, K , rj, rf are important for our understanding of the factor¬ 
ization hypothesis and of the origin of the nonfactorizable contributions. Namely, there is a 
huge discrepancy between the experimental results for some of the decays and the theoreti¬ 
cal predictions based on the factorization. Even the inclusion of calculated nonfactorizable 
contributions in some of B —> X cS K decays [2] has not shown satisfactory agreement with 
the experiment. Recently we have extracted the decay constants of charmonia states by 
LCSR and by the lattice calculations [3]. With the determined form factors of transitions 
B, B s —y rf 1 ' 1 in this paper it will be possible to analyze consistently nonleptonic decays to 
charmonia and to test the factorization hypothesis in such transitions. 

Decays B s —>• X c5 P are also useful to access CP violation in the B s sector and the 
phase of the B s — B s mixing, f3 s = arg (—Vt s V t *b/V cs V* b ) [4] and in the combination with 
the B —> X C qT)^ observables they can be also used for the determination of the rj — r/ 
mixing parameters [5, 6]. 

By using the huge amount of data it could be possible to make a thorough analysis and 
to extract the nonfactorizable contributions of nonleptonic decays from the data. The first 
ingredient for the analysis is certainly our knowledge of the B(D) —> P and B S (D S ) —>■ P 
form factors. These form factors have been calculated for years by using the QCD light- 
cone sum rule (LCSR) method [7] and on the lattice, step by step improving the precision 
of the results. The form factors for B(D ) —>• 7r, it and B S (D S ) —> ir,K are known now 
with quite a remarkable precision due to the consistent inclusion of corrections up to the 
twist-four a the LO and up to the twist-3 at the NLO [8—11]. 

With the recent update on the rj, rj' DAs where the SU(3) breaking effects are included 
consistently to the power-suppressed twist-four corrections [12], it is possible now to analyze 
B(D ) —> t/ 7) and B S (D S ) —> r /7) form factors to the same precision as for the B{D) —> 7r, K 
and B S {D S ) —> tt,K. But, ?/ and rf mesons exhibit some issues which makes them quite 
different form the pion. In the exact SU(3) flavor limit rj is a pure flavor-octet state, 
while rf is a pure flavor-singlet. Due to the existence of the axial 17(1) anomaly, i.e. the 
SU(3) breaking effects which are large and responsible for the heaviness of rf , there is a 
mixture between flavour-octet and flavour-singlet states usually described by the mixing 
matrix. In addition, the flavour-singlet states can mix with the two-gluon states producing 
the large gluonic admixture in rf mesons (which are primarily flavour-singlet states) and 
almost negligible ones in rj mesons. These gluonic contributions to the B(D ) —> t/7) and 
B S (D S ) —> rjW form factors enter at the NLO level which makes them quite nontrivial for 
calculation. The only existing calculation was done by Ball and Jones [13] for the f + form 
factor of the B —> rj ' > decay. 

We check those results, improve them by including the m (,) corrections to the both, 
the hard scattering amplitude and to the DA of 7/7) and consistently combine them inside 
the rj — rf mixing schemes with the ’standard’ quark contributions to predict B —> rj 7) 
but also D —> 7/7) and B S (D S ) —> t/ 7) transition form factor / + . In order to calculate 
consistently rare semileptonic B(D ) —>• 7/7) and B S (D S ) —> 7/7) decays such as, for example, 
B s —y t/ 7)/+/ _ and B s rf'^vv, it is necessary to calculate also other form factors, f° and 
f T (for definitions see (3.2, 3. 3, 3. 7)) of these decays which is for the first time done in this 
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paper. 


2 r) — rf mixing schemes and distribution amplitudes 
2.1 Mixing 

To analyze r\ and rf states, we have to deal with several definitions of matrix elements of 
the flavour-diagonal axial vector and pseudoscalar current: 

(P(p)|g7 M 75<?|0) = -~^=f q p p^, 2m q (P(p)\q~f 5 q\0) = ~^ h P’ 

(P(p)|s7^75s|0) = - ifpp M , 2m s (P(p)|s7 5 s|0) = -ih s P , (2.1) 

where q = u,d and the isospin limit is taken, rn q = \{m u + rrid)- There is also a U(1)a 
anomaly, 

(F(p)|gG7C 4 -'“'|0) = ap, (2.2) 

which is connected with derivatives of the currents through the equation of motion as 

Wrf'lbQ) = 2im q q'y 5 q - ^G A V G A ^ (2.3) 

and included in hp S as 

ap= 7^|m| =h , p _ f , m2p (2.4) 

In the exact SU(3) flavour-symmetry limit ap = 0. 

It is known that the SU(3) breaking corrections for rj and rf are large and that i] and 
rf mix since they are not pure a flavour-octet and a flavour-singlet states, respectively. 

The mixing of r/ and rf mesons is established in two mixing schemes: the singlet-octet 
(SO) and the quark-flavour (QF) scheme. Each of the schemes has some advantages and 
some disadvantages. 

In the SO scheme the mixing occurs among SU(3)f singlet 1 771 ) = l/\/3| uu + dd + ss) 
and octet 1 77 s) = l/\/6| uu + dd — 2 ss) components. By defining the coupling of the axial- 
currents to r] and rf mesons as 

(p)) = ify)Pn , (i = 1,8), (2.5) 

the decay constants of pure (hypothetical) singlet and octet states f % are related to the /*, 
via two-parameter mixing matrix 

(4 fv\ - ( COS08 ~ s ' m9 l \ (fs°\ 

V/J fry) \sin 08 COS 9, J { 0 f\) ' 1 ‘ ’ 

Since only singlet component mixes with the gluonic contributions, the renormalization 
scale dependence of parameters is diagonalized in the SO scheme and therefore is suitable 
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for the analysis of the gluon distribution amplitudes [14]. Moreover, /§ is scale independent 
and fi renormalizes multiplicatively: 


/!(/*) = /p(m o) 

/p(^) = /p(mo) (i + ^ Mm) - « s (mo)]1 , (2.7) 

where /to = 1 GeV, the scale at which the values of the mixing parameters are determined 

[15]- 

The simpler mixing scheme is QF scheme. There the basic components are \rj q ) = 
1/V2\uu + dd) and | r} s ) = |ss) states and the decay constants are defined as 

(p)) = ify)Pn , ('r = q,s ). ( 2 . 8 ) 

Their mixing with the decay constants of pure (hypothetical) non-strange and strange 
states, f q and f s respectively, is given by 

(fv fv\ = ( cos 6q -smd s \ (fq o\ (2q) 

\.fq> fq' j \sin e q cos6» s J yo f s J 

The main advantage of this scheme is that the mixing is not governed by the (large, 
10 — 20%) SU(3)f breaking effects as in the SO scheme, but by the OZI-rule violating 
contributions which have be proven to be small [15]. Therefore it is possible to parametrize 
the mixing just with one angle </> and the matrix U{4>) given as 


u(fa = 

which leads to the following expressions 


cos o — Sill ( 


( 2 . 10 ) 


V 


SSH'TT 

fn q) = fq COS (f> = (\/2 fq + /®) . 

= fq sin (j) = (V2f v , + 


= U(fa 


Vq 


fri &) = ~fs sin / 
fif 1 * = fs COS 4> 


The parameters have been determined by fits in [15] as 
f q = (1.07 ± 0.02)/, , f a = (1.34 ± 0.06)/, , 


( 2 . 11 ) 

= (7S W ~ V^/, 8 ) . 

= (75 {fi - ^ff) ■ 

= 39.3° ±1.0°, (2.12) 


and will be also used in this paper 1 . These values give for the parameters of the SO basis 
the following: 


f 8 = (1.26 ± 0.04)/, , f\ = (1.17 ± 0.03)/, , fa = -(21.2° ± 1.6°), fa = -(9.2° ± 1.7°), 

(2.13) 

1 There have been some recent discussions on the 27 — 17' mixing parameters and all of them are in the 

range of (f> above [6, 16 - 18 ]. 
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and the decay constants are connected as 





/ 0.1530 0.0243 
0.0595 0.1506 


GeV. 


(2.14) 


Due to the mixing of the flavour-singlet and gluonic components, in the QF scheme both 
r) q and ij s will get gluonic contributions and therefore also the physical ij and rf states. The 
flavour states in QF scheme and in the approximation above can be written as [14] 

W) = {^q( x ^F)\qq) + y/2/3il> g (x, fi F )\gg)) (2.15) 

> = 2 jiw + V^T n F )\gg)^ (2. 

where \qq) = (uu + dd)/yj 2 and ?/> q = 1/3(i/js + 2 - 0 i) and ij) s = 1 / 3 ( 2-08 + V’l)- 

By combing above information about the nature of rj and rf states one can expect 
that gluonic contributions \gg) will be larger for r/' mesons, which is confirmed by the final 
results. 

Until now there is no available QCD sum rule or lattice QCD calculations of B s to 77 ^ 
transition form factors /+’^. Since these transitions probe only the |ss) content, one can 
use the approximation in the quark flavour scheme 


fB s ri = - Sin 4>fBK > fB sT]' = cos 4>fsK , ( 2 -17) 

which neglects completely the gluonic contribution. The calculation presented in this paper 
will check for the SU(3)f breaking effects in the above relations. 


2.2 Distribution amplitudes 

The light-cone distribution amplitudes (DA), giving the momentum fraction distribution 
of valence quarks of rj and rj' are defined analogously to other meson light-cone DAs, by 
expanding the non-local operators on the light-cone in terms of increasing twist, but paying 
attention to the specific flavour structure of 77 W mesons. 

The twist 2 two-quark DAs ^ p °f F = r/ / ) mesons are defined as 

< 0 |*(*)C^75[*, -zM-z)\P(p)) = i(pz)fp T (u) , ( 2 . 18 ) 

Jo 

where as usual is the light-like vector and [z,—z] is the path-ordered gauge connection 
and u is the momentum fraction of a valence quark. In the SO basis one will have C\ = \j \/3 
and Cg = \g/\/2 ( As is the standard Gell-Mann matrix ), while in QF basis the constants 
are C q = (\/2Ci +Cg)/\/3 and C s = (Ci — \/2Cg)/y/3. The twist -2 two-quark DAs of 77 ^ are 
symmetric in their argument and therefore they can be expanded in terms of Gegenbauer 
polynomials as usually: 


W 2,P 


(u) = 6 u(l - u) 1 + a n t (f i ) c n 2 ( 2u - 1 )] (* = 1,8 ,q,s). (2.19) 


n=2,4,.. 
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Pi 

The coefficients a n ’ are the Gegenbauer moments of the quark DA. 

The gluonic twist-2 DA 0® P of P = rp' 1 mesons are defined by the following matrix 
element (for detailed discussion on the derivation of gluonic DA and its mixing with the 
quark states in mesons see for example [19]): 

n^{0\G^{z)[z,-zK(-z)\P(p)) = \^(pPf'p j( 1 <i Ue ‘< 2 “- 1 ><>”>4>» F (u ) . (2.20) 
It is antisymmetric and therefore 


02,p(“) = -<^ P ( 1 - u ), 

5/2 

and it is expanded in terms of C n ‘ Gegenbauer polynomials 


( 2 . 21 ) 


$j,p(«) = « 2 (! - u f ( b n 9 {v) C n-l( 2u - !) 


( 2 . 22 ) 


, 71=2,4,.. 


where the coefficients b B ' 9 are the Gegenbauer moments of the gluon DA and we take 
bn 9 = bn ' 9 and keep only the first term in the sum, n = 2. Although bn 9 and bn ’ 9 could 
differ, this approximation is justified since their values are subject of large uncertainties. 
In the calculation we use the following matrix element of the over two gluon fields 




1 zPp a C F tl 5 ab 



due^ pz) 


u( 1 — u ) 


(2.23) 


With the above normalization of the DA, the renormalization mixing of twist-2 quark 
and gluonic distribution amplitudes is given as 


d (af )A \ 
% \bf’ 9 ) 


a s (d) 

4ir 


100 _io\ 

9 81 ] 

-36 22 J 



(2.24) 


and it is numerically small. But, the mixing is important for p 2 = 0 case, since it verifies 
the collinear ’factorization formula’ for the form factors 

rl 


F(q 2 ,(p +q) 2 ) = du^T { H\u,q 2 ,(p + q) 2 ,pi R )'i/j n {u,p m ), (2.25) 

40 


and proves that the separation of the transition form factors in perturbatively calculable 
hard-scattering Tp part and a nonperturbative DA is essentially independent on the fac¬ 
torization scale /iiR [20] . This is an essential step of calculation which is going to be proved 
for each of the T-correlation function at the order of twist 2, see discussion in the next 
section. 

The explicit solutions of (2.24) can be find in [13] and in the appendix B of [12]. 

In the asymptotic case, when Q 2 = — q 2 —> oo the twist-2 quark and gluon DAs evolve 
to their asymptotic forms 


</>2,p(' u )|asym = 6tt(l - u) , 

^2,p(' U )|asym = 0. (2.26) 
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In that case, there is no gluonic contribution at the twist-2 level to the form factors, and the 
residual /im dependence in the twist-2 NLO quark contribution integrates with (f> l 2 p(u)| asym 
to zero, which again confirms the fj,jp independence of the complete result. 

To include SU(3) flavour-breaking corrections consistently we keep not only cor¬ 
rections and quark masses in the hard-scattering amplitudes, but also in the distribution 
amplitudes. Therefore we do not use the approximations in the twist-3 and twist-4 contri¬ 
butions employed in the literature where the following replacements are used in DAs: 


h 


ml 


U 


ml 


h 


mz 


2 m| - ml 

JS 


(2.27) 


for M —) 
scheme): 


2?7ig 2 m q 2 m s 2 m s 

rj q and M —> r] s decays respectively. Instead we are going to use (in the QF 


O / o o o 9 /— / O o 

—>• h q = fq{mt cos (j> + mti sin <f> ) — V2 / s (mi — ml) sin<^cos0. 


Uml 


fq 


h s = fs {mil cos 2 <fr + m q sin' 4>) - 7 = {mil — m 2 ) sin cj) cos ( 

v 2 


(2.28) 


Although the above quantities, especially h q , are weakly constrained due to the numerical 
cancellations, 


hg = 0.0015 ± 0.0040 GeV 3 , h s = 0.087 ± 0.006 GeV 3 


(2.29) 


we use them for the consistency of our calculation. Actually, we will see later that the 
approximation in (2.27) for h q is quite bad and causes somewhat large values of form 
factors of D, B -» 77 ^. 

Distribution amplitudes of higher twist are defined following [12] and [21]. Their 
parameter evolutions and definitions include now the anomaly contribution ap with the 
following expressions [22]: 


a v — ^^{fqmq 
a v' — ^2 {■'i rn v' 


h q ) COS < f > = 

— hq) Sliuj) 


2 2 
m, — m 


v' ir i ( r~ \ 

-sin 4> cos 4> y—fq sin (f> + V 2 f s cos (j)J , 

sin 4> cos 4> ^ f q cos cj) + V2 f s sin cj^j . 


m-, - m- 


V2 


(2.30) 


Therefore in [12] the normalizations of two-particle twist-3 DAs differ from those in 
[21], In [12] one can find a consistent treatment of m s corrections up to twist-4 and of 
anomalous contributions to DA and we take definitions and expressions given there. 

Then, 

2m T (0\f{z2n)i'y5r{zin)\P{p)} = / due~ l ^ UZl+UZ2 ^ pn ^ <^p P {u) , (2-31) 

Jo 

where r = q,s and 

2m r (0\f{z 2 n)a llu 'y 5 r{z 1 n)\P{p)) = ^ " 2 \ p^n u - p^n^) [ e~ l ^ uzi+uz2 ^ pn) cj)^{u). 

o Jo 

(2.32) 
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The normalization is then 


where 


and 



(«) 


du(t> 3 p a (u) = H^p ] 


Hp ] = m 2 pFP - op 


H^ = 


h^_ 

V2 


= h 


(s) 


,. M) 
M) _ Jp 

P y/2 ’ 



f ( s ) 

Jp 


(2.33) 


(2.34) 


(2.35) 


By calculating the mixing of twist-4 DAs, some approximations in the twist-3 DA are 
made in [ 12 ] when compared to the expressions in [ 21 ], to keep the same order of calculation 
in the conformal spin and the quark masses. 

The expressions for the two-particle twist-3 DAs used (contributions of higher confor¬ 
mal spin and 0(m 2 ) corrections are neglected; see also [21], Eqs.(3.25-26)) are 


= 6 «(! 


4>l s = h s + 60m s / 3s C , 2 /2 (2M 
-u) h s + 10m s f 3 s Co / 2 {2u - 


- 1 ), 

1 )" ■ 


(2.36) 


The three-particle quark-gluon-antiquark DA is defined as usual [21] 


(0\r(z)a flu 'y 5 gG a p(vz)r(-z)\P(p)) = if 3r [p a p ^3 ~ PaPugjjp ~ (a /?)) 


daida 2 da 3 5(l - aq — a 2 — a 3 )® 3 r (ai,a 2 , a 3 ), (2.37) 


d’3r(a) 


360a±a 2 a 3 |l + A 3r (ai - a 2 ) + u 3 r -(7a 3 


(2.38) 


There are two two-particle twist-4 DAs V’i 'p{ u )i ^ 4 T p( u ) and four three-particle twist-4 
DAs, ’kjp(a), T^p(a), <h^p(a), &^p(a). All details and subtleties in derivation of these 
improved twist-4 DAs with the corrected mass corrections and inclusion of the anomalous 
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contribution can be found in Appendix A of [12]. Here we just quote the expressions: 

</$(«) = d’p’» + 

,$‘ w (u) — — 5 2 p^ cy 2 (2u — 1) + 30 ( — — 10i.(l — u) + 35u 2 (l — u) 2 ) , 

/ D \ ' 


s\ 

| - i9«d - u) + - ») 2 +4? 


- - 54m(1 -u) + 225m 2 (1 - uf 


02 w = < 2 tw (<‘)+'"^ 2 “ ! ('‘), 

42 ,w (u) — -!]^3p(r)ii 2 (l - u) 2 + 21 5 2 p r ^ujfp [u(l - u)(2 4- 3u(l - u)) 

O 

+2(m 3 (10 — 15m + 6m 2 ) lnw + (m H (1 — m))] 

Ar) 

+20m r -yM(l — u) (12 — 63m(1 — m) + 14m 2 (1 — m) 2 ) , 

fp 


42 m “(») = 1.(1 - u) 


88 39 , . 2/i \2 

— + —m(1 - m) + 14m 2 (1 - m) 2 
15 5 


and 


+ 


'p(a) = 120oiO 2O3 («i — a2)4>ip 


(r) 


<I , 4p(a) = 120aia2«3 (f>o/p + (303 — l)</>^ 


(r) 


I(r) 


2,P 


^4p(q;) = -30ag(ai - 02 ) 




(r) 1 

4 Pi 


V^p + + ^( 5 a 3 - 3 )v 4 "p 

(a) = - 30 o| (1 - a 3 )V’o r p + (03(1 - 03) - 60102) V’lp 


where 


t >( r ) - 
9 i,p — 


21 

¥ 


+ ( 0 : 3(1 _ 03 ) - -(a 2 + a|)) v 4 ,p 


<5p r ^ip + ¥ m P f 1 —p4p 


45 


-‘-2P 


o6 (r) - -^' 2(r) 
v'o.p — s°p ’ 

7(r) _ ^lf2(r) (r) 

p 2,P ~ g°P W 4P ) 

¥ r ) - _I<¥ r ) 

wo,p — 3 °p 1 


¥ r ) ~ 4 

,P “ 4 


¥ r ) - 4 

^2,P - 4 


5p (r) 4J + ^nip ^1 - y 45^ + 4m, 
25p (r) a;S - y ^1 - y - 4m,, 


fh) 

J 3P 
' f (r) 
JP 

fix) 
J 3P 

fO) 

JP 


— C2pM(l — M) -yM(l — M) + 180 m 2 (1 — m) 2 ^ 

——c^p^ [m 3 (10 — 15m + 6m 2 ) InM + (m <->• (1 — m))] , (2.39) 


(2.40) 


(2.41) 
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The parameters which appear here are parametrization of various local matrix elements 
and their values are taken from [21] and listed in Appendix B. 

The above twist-4 expressions are valid for flavour-octet contributions where there is 
no mixing with the gluonic twist-4 DA. For the flavour-singlet case one has to take this 
mixing into account. In the approximation taken in [12] the twist-4 DAs are given by the 
replacement 

-> hP = m 2 P fP - a P (2.42) 


everywhere at the twist-4 level where the mass mn 2 P occurs. As it was discussed in [12] this 
substitution ensures for the given accuracy the consistent normalization of the twist-3 and 
twist-4 DA and ensures that the same mixing FKS scheme applies also for the higher-twist 
contributions. 

For the values of parameters involved we will use crude estimates in terms of the pion 
and kaon DA parameters derived from the sum rules [21], see Appendix: 


P.q P.s 7r 
a 2,4 — a 2,4 — a 2,4 i 


f3q — fz-n i 

hq - 0 , 

U 3q ~ W 3 7r , 

K4q^0, 

A 2 D) ~ A 2 
Op — o^ , 


/3s — f3K , 

A 3s — A 3 r , 

U 3 s ~ U 3K , 


K^g 


K4K , 
s; 2 


A 2 ( s ) ~ A 2 
Up — u K , 


while the corresponding rj, rj parameters will be given through the mixing as 


0 


/ Aq) A s )\ (f 


4“ /,<•>' 

A- 


- u <*> Vo l 


(2.43) 


3 LCSR for B, B s —>• 77 W and D, D s —> 77 W form factors 

For calculating the M —> rf 1 ' 1 form factors, where M = B, B s , D, D s , by using the LCSR 
method one considers a vacuum-to -77 or vacuum-to-r/ correlation functions of a weak current 
and an interpolating current with the quantum numbers of a meson M. For B —> rf'\ the 
form factors w and /^ (;) will be defined with the help of the correlator 

Fp(p,Q) = i J d*x e iq ' x {rj ( ' ,) (p)\T{u(x)r iM b(x),j B (0)} | 0 ) 

f F (q 2 i (p + q) 2 )Pn + F (q 2 , (p + q ) 2 )%, = 7a* 

(3T) 

^ F T (q 2 , (p + q ) 2 ) [p^q 2 - q^qp )], r M = -ia pv q v 

for two different b —>• u transition currents, where jp = m^n^u. Analogous formulas are 
going to be valid for D —> 77 W with the replacements b —> c and u — > d in the transition 
currents and jp -> jo = m c ar(§d. For form factors we consider the replacement 
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Figure 1. Diagrams corresponding to the leading-order terms in the hard-scattering amplitudes 
involving the two-particle (left) and three-particle (right) rf'1 DA’s shown by ovals. Solid, curly and 
wave lines represent quarks, gluons, and external currents, respectively. For B s —> 77 W transition, 
u is replaced by s. In the case of D —► 77 ^ transitions, u —» d and b —» c and correspondingly d is 
exchanged by s for D s — > 77M. 



Figure 2. Diagrams contributing to the quark hard-scattering amplitudes at 0(a s ). 


u —>• s in (3.1) and js s = m^bi^s interpolating current. Again, D s case is then obtained 
trivially by replacing 6 -quark with the c-quark. 

Since we want to explore also the SU(3) symmetry breaking, we will keep the rp 1 ' 
masses ( p 2 = rnr w ) in (3.1). The light quark masses will be systematically neglected, 
except when they occur in ratios in the distribution amplitudes. 

The method of the LCSR is very well know and we will here just briefly outline the 
procedure in order to properly define all ingredients necessary for calculating the form 
factors. For the large virtualities of the currents above, the correlation function is dom¬ 
inated by the distances x 2 = 0 near the light-cone, and factorizes to the convolution of 
the nonperturbative, universal part (the light cone distribution amplitude (DA)) and the 
perturbative, short-distance part, the hard scattering amplitude, as a sum of contributions 
of increasing twist. 

We calculate here contributions up to the twist-4 in the leading order, 0(oP s ), and up to 
the twist-3 in NLO, neglecting the three-particle contributions at this level. Schematically, 
the contributions are shown in Fig.l and Fig.2. Due to the specific properties of 7/ and rf 
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Figure 3. Diagrams contributing to the gluonic hard-scattering amplitudes at 0(a s ).The first 
diagram is IR divergent and its divergence will be absorbed by the evolution of the gluon DA. See 
text. 


mesons discussed above, there are additional gluonic diagrams contributing to M — > rj-' 1 
form factors shown in Fig. 3. These contribution has only been calculated for /+ form 
factor at twist-2 level in [13] and for = 0. Here we are going to calculate these 

contributions for other form factors w and by neglecting 0(a s m v (i)) effects in 

both DAs and the hard-scattering part. This approximation is justified having in mind 
that parameters of DA for the gluonic DA of p and rf are badly known, see the values of 
&2 ,s parameter below. 

By using hadronic dispersion relation in the virtuality (p + q) 2 of the current in the B 
channel, we can relate the correlation function (3.1) to the B —> rj W matrix elements, 


{■q ( '\v)\uAp.b\B(p + q)) = 2/+^ (q 2 ) Ptl + (<? 2 ) + f Br](l) (<? 2 )) q„ , (3.2) 


{p {,) (p)\ua^q u b\B(p + q)) = q 2 (2p^ + q^ - ( m% - m 2 (/) )q 


J m B + m v o) 


(3.3) 


and extract the form factors. In the literature it sometimes appears that the form factors 
are defined as above by divided by a factor \pl to match the transition form factors of 
rf, rf with those of a pion when there is no rj — rf mixing and in the limit of the conserved 
SU(3)-flavour symmetry [13]. 

Inserting hadronic states with the H-meson quantum numbers between the currents 
in (3.1), and isolating the ground-state H-meson contributions for all three invariant am¬ 
plitudes F(q 2 , (p + q) 2 ), F(q 2 , (p + q) 2 ) and F T (q 2 , (p + q) 2 ) and using (3.2) and (3.3) 
obtains: 


/Xw(^ 2 ) = 


2 m 2 R f B 


F 0 (q 2 , M 2 , s*) + ^ ( F^q- 2 , M 2 , s$) + F? 9 ’ + (q 2 , M 2 , s*)) 


e m 2 B /M 2 


^BnO ( q ) + -W) ) m\f B 

T , 2 \ _ (™B + m^)e m s/ M2 
2 m 2 f B 


F 0 (q 2 , M 2 , s$) + ~^~Fi{q 2 , M 2 , s$) 


(3.4) 

(3.5) 


+ 


B- 

a s C B 

4 - 7 T 


Fq {q 2 1 M 2 , Sq) 

(f? ( q 2 , M 2 , s$) + F^ T ( q 2 , M 2 , sj?j) 


(3.6) 
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The scalar B —> 77 ^ form factor is then a combination of the vector form factor (3.4) and 
the form factor from (3.5), 


«(<Z 2 ) + 


m B — m 




M 


(3.7) 


and is only present in the semileptonic B^, decays when the lepton mass is 

not neglected and the rare B S ,D S —> decays. In above, F o(i) and F 0 (!) represent 

the LO (NLO) contributions and fs = (Bd\mbbi-f5cl\0)/rn“^ is the 13-meson decay constant. 
pWi ■.{++) are } eac [ing order twist-2 two-gluon contributions calculated explicitly in the pa¬ 
per. At the leading twist-2 level there is no gluonic contribution in (3.5). However, note 
from (3.7) that this does not mean that twist-2 two-gluon contributions will not appear in 
the scalar w form factors (3.7). 

As usual, the quark-hadron duality is used to approximate heavier state contribution 
by introducing the effective threshold parameter and the ground state contribution of B 
meson is enhanced by the Borel-transformation in the variable (p + q ) 2 —> M 2 . Completely 
analogous relations are valid for B s —> rp'> form factors, with the replacement u —> s in 
(3.2) and (3.3) and by replacing ms by mB sl Jb by fs B , as well as M 2 by M 2 and Sq by 
Sq s in (3.4 - 3.6). In addition, in the derivation of above expressions for B Sl one has to take 
into account that (B s \bir/ 5 s\0)= fB a /(nib + m s ). The same is valid for D,D S form 
factors with the replacement m^ —> rn c and the appropriate exchanges described before. 

The calculation will be performed in MS scheme. The B, B s and D, D s decay constants 
will be calculated in the MS scheme using the sum rule expressions from [23] with O(o s , m 2 ) 
accuracy. In that way we achieve the consistency of the calculation and the cancellation of 
uncertainties in the sum rule parameters. 

Each form factor can be written in the form of the dispersion relation: 


*o 

F(q 2 , Mf s) ,sft 1 ) = i J dse~ s/M M lm s F(q 2 , s ), (3.8) 

ml 

where now s = (p+q) 2 ■ 

The leading order parts of the LCSR for /+ )?(;) , + f M ^, ) and form factors 

are given in Appendix A. 

Up to now, SU(3)-violating effects for f D (,), f B (/ ) form factors were not system- 
atically studied, since the effects of inclusion of m 2 w effects complicate the calculation, 
especially at NLO in the hard-scattering amplitudes. However, while the complete SU(3)- 
symmetry breaking corrections in rf'^ DAs of twist-3 and twist-4 are now known [12], it 
is worth to have a consistent picture of all SU(3)-breaking corrections and we will include 
complete SU(3)-breaking effects in both DAs, as well as in the hard scattering amplitudes 
at LO. At NLO in the hard-scattering amplitudes, for the cases when the mass of a light 
quark cannot be neglected, as for m s , the inclusion of m s and w- 2 (/) effects complicate 
the calculation. As already known from the analysis of B^ —> K from factors done in 
[11], inclusion of quark mass effects leads to the mixing between different twists and the 
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fully consistent calculation with m s included in the quark propagators is not possible, see 
discussion in [11], However, here we have rf as a finite-state particle which mass is much 
larger than m s and therefore, in the NLO quark and gluonic amplitudes we set m s = 0 
and p 2 = m 2 , / 0. 

Each form factor can be expressed as 


and 


' f +,0,T 
J Br) 


\t = U(<t>) 


I Br/' 


/ r+,0 ,TN 
J Bril 

jBr) s j 


(3.9) 


r+,o,T _ Aqq) +,o,^ A 99 ) i*+,o,T _ A 99 ) 

J BtjQ Br ' J BrjQ ’ J Brj s J Brj s 


(3.10) 


and explicitly 


f +,0,T 
J Brj 


Aq) 

h_^ + F f q )+f^Ff MT 


Aq) 

f^ T =^W + Fr)+f^ 9 ^\ 

= fi s) {^ s + Ff s ) + flFfW , 

f+,0,T r(s) / rpss 1 pss\ 1 rl rpgg,+,0,T 

JBs-q' - V cO + *1 ) + Jr,'* 1 ! 


(3.11) 


where F^ q and Fq s ( F qq and F ( s ) are LO (NLO) contributions from quark hard-scattering 
amplitudes for each of the form factors and F gg is the NLO gluonic contribution propor¬ 
tional to the singlet-flavour decay constants 


f 1 = — 
/r? v/3 

fri ' V3 


2 COS (j)f q 



2 sin (j)f q + cos 4>f s 


(3.12) 


The decay constants are given in ( 2 . 12 ). Analogous expressions are valid for D^ —> 77 W 
decays. 

Obviously, for B,D —> 77 W transitions the main contribution conies from p q me¬ 
son states and p s contributes only through suppressed gluonic contributions, while for 
B S ,D S —> rf'' transitions the leading rj s meson state contribution will receive, through the 
gluonic diagrams, a small mixture with p q state. Also, implicitly there will be mixing with 
among twist-2 quark and gluonic distribution amplitudes Eq.(2.24), which will bring hlf ^ 3 
dependence in the twist-2 quark LO {Fq Q and Fq s ) and NLO contributions (F qq and Ff s ) 
and a g ^ dependence to the gluonic contributions F gg . 

Since rj and rf are mixtures of the | qq), |ss) in the calculation of the quark contributions 
we will use (with appropriate substitutions) our NLO results for the hard-scattering part 
for B —> 7 r [10] and B^ —> K form factors [11] with the p 2 effects included at the LO 
(up to twist-4) and NLO level (up to twist-3) and will imply recently derived DAs of p 
and rf with the SU(3)-breaking effects and the axial anomaly contributions included up 
to twist-4. The gluonic contributions, which are already NLO effect, will be calculated for 
V 2 = m 2 ^ = 0. 
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4 LCSR for gluonic contributions to the form factors and consistent 
treatment of the IR divergences appearing 

The gluonic contributions at the 0(a s ) to the B(D ) —> 77 W and B S (D S ) —> rj W form factors 

1 0 J 1 

come from the diagrams in Fig.3. The results for the form factors are presented 

in subsection 4.1. They are added to the quark contributions (3.11) to get the complete 
result at the order 0 (a s ). 

The first diagram is Fig.3 is IR (collinear) divergent. This divergence has to disappear 
for the general collinear factorization formula used here 


F{q 2 ,{p + q) 2 ) = ] (u,q 2 ,(p + q) 2 ,ii F ) ® & n , P (u, fi F ) . 


= [ du. 

Jo 


(4.1) 


be valid. The scale is the factorization scale. At the twist n = 2 level, as already 

mentioned, there will be mixing of quark and gluonic contributions and the hard-scattering 
( 2 ) 

(perturbative part) T H ; and the distribution amplitude 4>2 can be represented as 


rp( 2 ) _ / Tqg 
\ T 99j 


4*2, P = 


b l,P 

^2 ,P> 


(4.2) 


In order to consistently treat this mixing we have to examine the evolution of the DAs, 
at the same 0(a s ) as the calculation of the perturbative part Tu. Due to the mixing the 
standard Brodsky-Lepage (BL) evolution equation [24] 

d 


Hf w— ip(u, [if] = V(u, v, /. l f ) ( 8 ) hf) , 
OHF 


(4.3) 


will be a matrix equation now, where V(u, v,fip) is the perturbatively calculable evolution 
kernel 

«I(a»f). 


V(u,v,hf) = —j^Vi (u,v) + 
with the LO kernel of our interest of the form 


(47r) 


-V 2 (u,v) 


(4.4) 


«(»,») = \l” 

V v gq v gg > 


(4.5) 


and V tJ are well-know evolution kernels [14, 19, 25] which we cite here for convenience: 


u 


V qq (u,v) = 2 C F < - 


1 + 


v — u 


Q(v — u) + 


u 


1 -u 

1 — V , 


Vqg(u, V ) = -2^UfC F { ^@(u -u) - 


u 


1 — u 


V gq (u, v ) = 2 y 'rifCp | —Q(v - u) - 

Q(v — u ) 


I 

Vgg(u,V ) = 2JV C j - 


v — u 


+ 


U -> 1 - r 

w —>• 1 — u\ 
v —> 1 — u y 

2 «- 1 


u 


0 (u — u) 


+ 


u 


1 — u 
1 -V 


+ /?o£(u _ 77) . 

(4.6) 
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These evolution kernels are exactly those which govern the renormalization of the DAs 

<f>(rx) = Z^ ren {u, v, lift) <g> ${v, fift) . (4.7) 

The connection between Z and the evolution kernel V is given as 

V(Tr) = ~ Z faren^Ft) {^R Z <t>,r e n{fijl)^ , Z cf>,ren(^R) = H ~Vi(u,v) H . 

(D = 4 — 2e). On the other hand, by calculating the hard-scattering part Th, owing to 
the fact that final-state quarks are taken to be massless and on-shell (for the case p 2 = 0), 
the amplitude contains collinear singularities. Since Th is a finite quantity by definition, 
collinear singularities have to be subtracted. Therefore, T factorizes as 

T(u,Q 2 ) = Th(v,Q 2 , up) ® z T,coi(v,u] pf) , (4.8) 

with collinear singularities being subtracted at the scale p f and absorbed into the constant 
z T,coi ■ As usual The UV singularities are removed by the renormalization of the fields and 
by the coupling-constant renormalization at the (renormalization) scale p,ft. Now, in order 
that the factorization formula is valid, the following has to be satisfied 

z T,coi(u,v; up) <g> Z<f, >ren (v,u; fj. F ) = $( v - w ) • (4.9) 

The divergences of T(u, Q 2 ) and 4>(tt) in (4.1) then cancel and at the end we are left with 
the finite perturbative expressions for all form transition factors 

F{q 2 , (P + q ) 2 ) = T h (u, Q 2 , hf) <g> 4>(u, hf) . (4.10) 

It is worth pointing out that the scale [ip representing the boundary between the low- and 
high-energy parts in (4.1) plays the role of the separation scale for collinear singularities in 
T(u, Q 2 ), on the one hand, and of the renormalization scale for UV singularities appearing 
in the perturbatively calculable part of the distribution amplitude 4>(rt), on the other hand. 
The general discussion and all details of the proof of the cancellation of the factorization 
scale dependence in the collinear factorization formula (4.1) at all orders of calculation can 
be found in [20, 26]. 

In our case of calculating the heavy-to-light transition form factors / +,0,T we face the 
following situation. The hard-scattering, perturbatively calculable pieces coming from the 
diagrams from Fig.2 have UV and infra-red singularities at 0(a s ). We have already proven 
in [10, 11, 27] for B —> ir and B ( s ) —> K form factors that the IR divergences of the quark 
contributions at twist 2 level cancel exactly with those coming from the evolution kernel 
V qq . Here, due to the mixing with the twist 2 gluonic contributions, the convolution of V qg 
of the Th LO will exactly cancel the IR divergence in the first gluonic diagram in Fig.3. At 
the twist 3 level of 0(a s ) the IR divergences of quark diagrams mutually cancel, as shown 
before in [10, 11]. This gives the final proof of the collinear factorization formula at the 
given order for the heavy-to-light M —> p' transition form factors. 
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4.1 Explicit results for the leading two-gluon contributions to the / + and f 1 
form factors in B,B S —> r/T and D,D S —> rf ' 1 transitions 

In the calculation of the gluonic contributions to the form factors we have faced the problem 
of the consistent treatment of the 75 in the dimensional regularization. Leading order for 
the gluonic amplitude is given by one-loop Feynman diagrams in Fig. 3. and we have to 
deal with IR divergence which is a consequence of having massless quarks propagating 
through the loops. In the calculation of the gluonic contributions to the form factors it 
appears a Levi-Civita tensor in the projector of the twist-2 two-gluon DA (2.23) and a 
single 75 matrix in the trace which are both quantities with well-defined properties only 
in D = 4 space-time dimensions. Generalization of these quantities in D dimensions is 
problematic and different approaches to avoid resulting ambiguities can be found in the 
literature. Moreover, in our case there is no gluonic contributions which appear at LO on 
a s that would greatly help in resolving the 75 problem at NLO level. The problem was not 
addressed in the paper where the gluonic amplitude was evaluated for the first time [13] 
and it is not clear how they resolved the ambiguities. 

In the case of the interest it is possible to completely avoid 75 problem and all connected 
complications since the IR divergence is direct consequence of the massless quark lines and 
putting a small mass m in massless quark propagators regularizes (removes) the divergence. 
As a consequence, we are not forced to use dimensional regularization and calculation can 
be performed in four dimensions without any problem. Note that putting mass in quark 
propagators doesn’t spoil any of properties and symmetries of the amplitude contrary to the 
case when, so called, mass regularization is used on gluon propagators. At the very end of 
calculation it is necessary to expand final result around zero for the small introduced quark 
mass m. The IR divergence will now reappear as ln(m 2 ) term and it is straightforward to 
connect it with l/(D-4) term in the framework of dimensional regularization. 

The obtained expressions are as follows: 


F? 9 ’\q 2 ,M 2 


s 0 ) 


_ fi L 

— J r ,(')°2 fj 


: 

J m z 


exp ""/" 2 f 99 ’ 


(», q 2 


(4.11) 


where the gluonic contribution to / + form factors is 


fgg,+ = 20m 


2 (« - m 2 ) 


27^3(s - q 2 ) 5 


3 (m 2 — q 2 ){hm A — bm 2 (q 2 + s) + q A + 3 q 2 s + s 2 ) 


2 In 


-lnlA, 


s — m 


m* J \ m 


- (37m 6 - m 4 (56 q 2 + 55s) + m 2 (l 8 g 4 + 76 q 2 s + 17s 2 ) + 3 q 6 - 27q A s - 11 q 2 s 2 - 2s 3 ) ^ 

(4.12) 
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and the corresponding contribution to f T form factors has the following form 

.2 


f 99 ' 1 = 5m- 


s — m 


12 q z 


m‘ 


s — m 


, ^ —9 (l 4 + 3g 2 s + s 2 ) In 

27\/3(s - q 2 ) 5 V s-m 2V ’ 

+ 6 (m 2 — q 2 ) (5m 4 — 5 m 2 {q 2 + s) + q l + 3 q 2 s + s 2 ) ^2 In ^ 

— (59m 6 — m 4 (72g 2 + 85s) + m 2 s(84g 2 + 23s) + 3 ( 6 < 7 6 + 6 q 4 s + 6 q 2 s 2 — s 3 )) 


- In ( ~~2 
m z 


(4.13) 


with m = m C) b- 

With respect to the fact there is no LO 0(a 6 ) twist-2 gluon contributions and following 
the discussions at the beginning of Sec.4, obviously there is no gluonic contributions to 
/ + + f~ form factors at this order of calculation. 

The result for the gluonic 0(a s ) contribution to f + form factors was first given in 
the appendix of [13]. Our result (4.12) does not completely agree with the one presented 
there. While we agree in the part being proportional to the logarithmic terms, there is a 
disagreement between the coefficients in the second line of (4.12) and the expression (A.l) 
from [13]. Since those terms are exactly those which change with the different treatment 
of 75 , and the authors of [13] have not placed any comment how they have resolved the 
75 ambiguities in the calculation of f"' + , we assume that the difference comes from the 
improper treatment of the 75 in [13]. 

The result for the gluonic 0(a s ) contribution to f T , Eq.(4.13) is a new result. 

5 Predictions for B, B s — > 77 ^ and D, D s — > rj^ form factors (/ + , f° and f T ) 

the form factors 

The prediction for B,B S — > 77 W and D,D S — > 77 W form factors (/ + ,/° and f T ) the form 
factors will be given in the MS scheme by using the input parameters listed in Appendix 
B. 

From expressions (3.5,3. 6 ,3.7) we see that we need the heavy-meson decay constants 
of B( s ) and in the calculation. As usually done, to achieve partial cancellation of 
the uncertainties in the calculation the two-point QCD sum rules for the decay constants 
/S)/b 0 an d fDlfDs is used in the same scheme, with 0(a s ,m 2 ) corrections included [23]. 
We have used the same level of accuracy as in the calculation of the form factors, i.e 
0(a s ) in both, the perturbative and nonperturbative (quark condensate) part and in the 
determination of the sum rules parameters s 0 <s) and M / s \ have used the usual consistency 
conditions in the sum rule calculations. 

The resulting predictions for /m, together with the fitted sum rule parameters for each 
of the mesons are given in the Appendix B, Tables 2-4. Here we quote the calculated values 
from Table 4: 


f D = 191 ± 9 MeV, f Ds = 219 ± 7 MeV, 

f B = 215 ± 7 MeV, f Bs = 246 ± 8 MeV, (5.1) 
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where the quoted error intervals are coming from the variation of and m\,j only since 
other uncertainties are canceled in ratios in Eqs.(3.5,3.6,3.7). By comparing our results 
with the previous LCSR results and the most recent determinations from [28], where in 
the perturbative part the higher order corrections were included, we see good agreement. 
The results are also within uncertainties of the lattice QCD calculations of the same decay 
constants [29], 

For the fo and fo s the experiment gives somewhat lower values [1], 
f D = 204.6 ± 5.0 MeV, f Ds = 257.5 ± 4.6 MeV , 

but still consistent within the complete LCSR error results [28]. 

The renormalization scale is given by the expression Hb^s) = — m b an d s ^ m ” 

ilarly for —> rf' 1 transitions. Therefore, for the renormalization scale we use /i = 3 

GeV, for the form factors and fx s = 3.4 GeV for and for /ip = 1.4 GeV and 

[iY) s = 1-5 GeV. As usual, we will check the sensitivity of the results on the variation of 
above scales and will include it in the error estimation. 

The method of extraction of the Borel parameters M and the effective thresholds 
So for /+JJ form factors is the same as described in [10]. It relies on the requirement 
that the derivative over —1/M 2 of the expression of the complete LCSR for a particular 
form factor, which gives heavy-meson masses rri 2 M , does not deviate more than 0.5 — 2.5% 
from the experimental values for those masses. Additional requirements such as that the 
subleading twist-4 terms in the LO, are small, less than 10% of the LO twist-2 term, 
that the NLO corrections of twist-2 and twist-3 parts are not exceeding 30% of their LO 
counterparts, and that the subtracted continuum remains small, are also satisfied. These 
demands provide us the central values for the LCSR parameters listed in Tab. 5. 

The estimated form factors for —> r /W are as follows: 

/£„(0) = 0.168lg;{$ = 0.168 ± 0.003 (b™) ± 0.002 (s 0 ,M) ±g;g& (mix) ±8:881 (rest), 

/X'(0) = 0 - 130 ±ao! = 0.130 ± 0.020 (bf 9 ) ± 0.002( So ,M) ±g;gg° (mix) ±8;88l (rest), 
1/^ (0)| = 0.212±g;gi| = 0.212 ± 0.003 {b v 2 ’ 9 ) ± 0.003(s 0 ,M) ± 0.012(mix) ±g;gg§ (rest), 

f+ sV ,( 0) = 0.252+°;° 2 ^ = 0.252 ± 0.019 ( b v 2 ’ 9 ) ± 0.004(s 0 ,M) ± 0.005(mix) ±g;gj£ (rest), 

(5.2) 

/*,(0) = 0.173lg;g^ = 0.173 ± 0.002 ( b± 0.003(s 0 , M) (mix) ±g;gg£ (rest), 
IbA 0) = 0.14llg;ggg = 0.141 ± 0.015 (bf 9 ) ± 0.002(s 0 , M) ±g;gl (mix) ±g;ggg (rest), 
l/L,(0)| = 0.225±g;gl| = 0.225 ± 0.002 (, b± 0.004( So , M) ±g;glf (mix) ±g;gjl (rest), 
fB sV >( 0) = 0.280±g;g?| = 0.280 ± 0.014 (bf 9 ) ± 0.004( So , M) ±8:88? (mix) ±g;8oi (rest), 

(5.3) 
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and for —> rj^: 


/d„( 0 ) = 

- 0 4 29 +0,165 - 

- u.+zy_ 0 141 - 

= 0.429 

± 

0.009 

M 9 ) 

10.004 

^O.OOl 

(so, M) ±0+41 

(mix) 

. 1 O .013 
±0.008 

(rest), 

= 

- 0 292 +0 ' 113 - 

- u.zyz_ 0104 - 

= 0.292 

± 

0.045 

{bt«\ 

V _i_0.009 
1 ±0.007 

(so, M) ±^ 9 ? 

(mix 

\ _i_0.015 
) ±0.011 

(rest), 

l/&„( 0 )l = 

- 0 495 + 0 ' 030 - 

- u.+yo_ 0 029 - 

= 0.495 

± 

0.007 

M 9 ) 

_i_0.004 

±0.002 

(so,M) ±8;81I 

(mix) 

. _l_0.016 

±0.009 

(rest), 

fit 77 ' ( 0 ) = 

- 0 f;^7 + 0 - 048 - 

- u.oo/_ 0 045 - 

= 0.557 

± 

0.041 

ibf 9 \ 

V 1 0.011 
> ±0.008 

(so, m) ± 8 : 80 ° 

(mix 

\ _i_ 0.022 

) ±0.014 

(rest), 











(5.4) 

/d»?( 0 ) = 

- 0 44 c ) +0 - 115 - 

- U.4oO _ 0 1Q7 - 

= 0.435 

± 

0.008 

W 9 ) 

_i_0.005 

±0.003 

(so, M) ± 0+06 

(mix) 

,0.177 

±0.151 

(rest) , 

bs 

"3 

- n QQ 7 + 0.118 _ 

- u.oo /_ ai47 - 

= 0.337 

± 

0.055 

{bi 9 \ 

V _l_0.013 
1 ±0.051 

(so, AO ± 8 : 1 S 

(mix 

\ _i_0.080 

) ±0.077 

(rest) , 

\fD s i)(®)\ ~ 

- n 441 + 0 - 091 - 

- u.++i _ 0 087 - 

= 0.441 

± 

0.007 

(bT) 

_i_0.052 

±0.005 

(so, M) ±° : °I? 

(mix) 

. _i_0.068 
±0.082 

(rest) , 

/d s 7j'(0) = 

- 0 655 +0 ' 072 - 

- o.ooo _ 0 065 - 

= 0.655 

± 

0.050 

(bt 9 : 

V _l_0.015 
1 ±0.014 

(so, M) ±o;oio 

(mix 

\ _i_0.034 
) ±0.026 

(rest) . 


(5.5) 


These results are predictions given with b 9 = 0 and then varied within the interval 
Abf 9 = ±20, which dependence is explicitly displayed in the errors. The errors are 
compilation of the variation of parameters added in quadratures. In the errors we explicitly 
stress SR parameter dependence (so, M), rj — rf mixing parameter dependence (mix) and 
dependences coming from the variation of the rest of parameters (rest={/i, m C} b, 02, 04}). 

The errors of the results are much larger for the transitions B, D -+ where B, D — > 
rjg dominates then for B S ,D S —> rj s decays since the error in the parameter h q (2.29) is 
huge, of 0 ( 200 %) depending not on r] — r / mixing parameters but exhibiting a numerical 
cancellation among terms. If one would use approximation (2.27) applied in [30] instead, 
the (rest)-errors would be almost an order of magnitude lower and the mean values would 
be somewhat larger for those decays, which we assume is the main reason, apart from the 
rest of SU(3)f approximations used there, of the discrepancies with some of the results 
presented in [30]. We see that the dominant errors in M —> rj form factors is coming 
from the variation of ' 9 and it amounts to about 15%, while in M -+ 77 decays come 
to 2%. Our findings for calculated B —> rf'^ form factors agree very well with those from 
[13, 31-33]. 

Their ^-dependence of the form factors and their ratios is shown in Fig.4-9. 

From Fig. 5 and Fig . 8 we see that the gluonic corrections are much larger for R( s ), D^ -+ 
rj decays then for M —> r/, as expected. Also the gluonic corrections are larger in D ( s ) 
decays. It is obvious that even in ratios of form factors the gluonic contributions give 
main error and that it would be difficult to constrain 62 , unless all M —> rj') semileptonic 
transitions are measured, Fig. 6 and 9. 

We can now investigate SU(3 )f approximations from (2.17). By using the obtained 
results and the result for Jbk from [ 11 ] we obtain 


I ±>+ I calc 
\JB s r ] \ 

/*+,calc 
J B S 77' 


0-212^3 vs 

0.252" vs 


l/s s r? 

f+> a 

J B sV ' 


approx I 


I i-)-,approx| 

\JE 


0.225 


+0.032 
-0.026 1 


0.278 


+0.038 
-0.031 ' 


(5.6) 
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Figure 4. Form factors for D ( s ) —> 77 W decays and their ratios. Solid lines represent f+ (/) 
/orm factors, dashed-dotted line )j)(/) and dashed line f^ jj(;) /orm factors. 

We can note that the approximation works quite well although somewhat better for M —> r/ 
decays than for M —> rf transitions. 

There exists LCSR calculations of f^ sT) form factor [34, 35]. In these papers the /^ , 
form factor is then obtained by using the relation 

= cot cj). (5.7) 

JD s ri 
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Figure 5. Gluonic dependence of f+ form factors. Shaded areas show change of the form 

factors under the variation of ' 9 = 0 ± 20. Solid line denotes the result for ’ 9 = 0, dashed- 
dotted for f >2 ’ 9 = 20 and dashed line for b ’ 9 = — 20 . 


While their predictions for agree with ours, the use of the above approximative relation 
which neglects the gluonic contributions gives somewhat larger , form factor then the 
one obtained here, (5.2, 5. 4). 

There exist also recent lattice results on D s —> rp ' 1 form factors [36]. These transitions 
at the lattice are challenging due to the presence of disconnected quark-line contributions 
and in [36] only the scalar f° D w form factors are calculated, which at q 2 = 0 are equal to 
the f + . By comparing the results one can see that the lattice predictions give /^ , < fo, : rp 
which is just opposite in LCSR for all M s —t transitions. The tendency in 

LCSR is established for non-strange meson decays, see results in (5. 2 , 5. 4). 


6 Phenomenological applications 

In this section we comment on some phenomenological results for semileptonic D—> rp 1 ’ 
and £>(<.) —>• 77 W decays which include the calculated from factors. To be able to calculate 
the branching ratio we need the form factor extracted in whole accessible kinematical 
regions. For decays the LCSR are applicable only in the region q 2 <C m 2 and for 
the region is 0 < q 2 < 12 GeV. 

The are many parametrization for calculating the shape of form factors at q 2 7 ^ 0. All 
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Figure 6. Gluonic dependence of ratios of D ( s ) —» ?/ , ) form factor ratios. Shaded areas show 


(0 


change of the form factors under the variation of ’ 9 = 0 ± 20. Solid line denotes the result for 


„(') 




(') 


&2 ’ 9 = 0, dashed-dotted for b'l ’ 9 = 20 and dashed line for b 2 ’ 9 = —20. 


of them work equally well and therefore we decided to use the most simplest one [37]: 


fM v 0)(r )|fit /«>) (i _ q*/m 2 H ,)(l - a+’ T q 2 /m 2 Ht ) 
f° Mv w(« 2 )|fit = (°)(i _ a o qZJmjj*) 


( 6 . 1 ) 
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q 2 ( GeV 2 ) 


Figure 7. Form factors for B^ —> 77 W decays and their ratios. Solid lines represent f+ form 
factors, dashed-dotted line |j(;) and dashed line fg ^ v w form factors. 


where the extrapolation of the form factors is performed just by fitting one parameter a 1 
for each of the decays and using the appropriate vector meson resonances rn* H , Tab. 6 , while 
the normalization is given by the form factors at q 2 = 0. The fitted parameters a 1 for 
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Figure 8. Gluonic dependence of w form factors. Shaded areas show change of the form 

factors under the variation of ' 9 = 0 ± 20. Solid line denotes the result for b% ’ 9 = 0, dashed- 
dotted for J >2 9 = 20 and dashed line for ’ 9 = — 20. 


form factors are 


a+ v = 0.165 ± 0.006, = 0.19 ± 0.05 

a+ sri = 0.198 ± 0.005 , a + , = 0.20 ± 0.03 (6.2) 

while for B^ are as follows: 

a Bri = 0-462 ± 0.002 , a° Bri = 1.00 ±0.01, a T Br] = 0.494 ± 0.005 

a Br) , = 0-45 ± 0.02 , oP Blf = 1.00 ±0.09, = 0.47 ± 0.04 

a B s ri = 0-505 ± 0.003 , a° Bgr) = 1.08 ±0.01, = 0.193 ± 0.002 

a B s rj' = 0-433 ± 0.003 , a° BsV , = 1.09 ±0.05, a^ slf = 0.51 ± 0.02 . (6.3) 


The semileptonic —> rf^eu e and B —> rf^ev e decay rates are calculated by 


( 0 , 




L H Jmi 


dq 2 \ 3/2 (q 2 )\f+ v(l) (q 2 )\ 2 , 


(6.4) 


where H = D,D Sl B and A (q 2 ) = (rn 2 H + rn 2 (l) - q 2 ) 2 - 4 m 2 H m 2 Q) and V Qq = V cd . V cs , V ub 
depending if D + , Df or B + meson is decaying, respectively. Values for the CKM matrix 
elements are taken from [1]: V c d = 0.225, V cs = 0.973, V ub = 0.0035, Vt s = 0.0405 ,V tb = 
0.999. (for V ub we used newly determined average value also from [1]). 
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Figure 9. Gluonic dependence of ratios of Br s ) — ► r/d form factor ratios. Shaded areas show 


(') 


change of the form factors under the variation of b 2 ’ 9 = 0 ± 20. So/zd /me denotes the result for 


„(0 


„(o 


(') 


&2 ’ 9 = 0, dashed-dotted for ’ 9 = 20 and dashed line for b 2 ’ 9 = —20. 


For the rare —?• (yv) decays we use the effective Standard Model hamiltonian 
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for b —> sl + l (w) transitions [38] and calculate decay rates as [39] 


G 2 \V fh V*\ 2 a 2 


V(B. -> ,<'»( + r) = , 

01Z7T Uio J4m 


dq 2 \ 1 / 2 {q 2 )\l- 


4 mf 1 


g 2 3g 2 * 


^ri r ) (.Q 


(6.5) 


where 


V>( ?2 ) = - m 5w) 2 | c 'io(^)/s 3 „(/)(9 2 )| 2 + (r-4m z 2 )A(g 2 )|Cio(^)/+ (/) (g 2 )p 


+(g 2 + 2mf)A(g 2 )|C 9 (/r)/+ riW (g 2 ) - 2 


TUb 4~ Tn s ^ t \ iT r 2\\2 


B ^ {,) ' J ' m Ss + m (/) 




and 


r(B s -> ?? (/ W) = 3 J Q S ’ ?() d? 2 A 3/2 (9 2 )|/+ >j(#) (? 2 )| 2 J (6.6) 


where Cl = Ci?/v / 2a/(27rsin 2 0\y)VtbV t *r]xX(xt) [39]. For the Wilson coefficients we use 
the following values 


C 7 = -0.3031, C 9 = 4.1696 , C w = -4.4641, C L = 2.74 • 10“ 9 . (6.7) 

Our predicted branching ratios for various M —> rj^ decays are given in Tab.l. By 
comparing with the existing calculations performed in different models [39, 44-46] we agree 
quite well, expect that we predict somewhat larger branching ratios for B s —> ?/ / V + r _ 
decays. 

Because of the larger errors in B,D —>• rf'^ decays, M s —y would be better for 

extraction of the unknown v</ parameter, but measurements of M s decays still have 
to achieve sufficient precision, in particular Br(B s —> ri'l + l~ ) and Br(B s — > r]l + l~ ) are 
challenging with the branching ratio of 0(1CC' — ICC 8 ) but they could be measured at 
future SuperB and SuperKEK experiments. 


7 Summary 

We have investigated B,B S — > rj M and D,D S —> rj M form factors (/ + ,/° and f T ) by 
including corrections in the leading (up to the twist-four) and next-to-leading order 

(up to the twist-three) in QCD, as well as gluonic contributions to the form factors at 
the leading twist in the framework of the QCD light-cone sum rules and have also taken 
SU(3)-flavour breaking corrections and the axial anomaly contributions to the distribution 
amplitudes consistently into account. The two-gluon twist-2 contributions are calculated 
for all f + , f° and f T form factors. 

We have given the values and shapes at q 2 ^ 0 of all calculated form factors and have 
shown predicted ratios for some semileptonic B, B s —> rj^ and D, D s — > rf'^ decay modes. 

With the determined form factors of transitions B,B S —>• rf'> it will be possible to ana¬ 
lyze consistently nonleptonic decays to charmonia and to test the factorization hypothesis 
in such transitions which we be a subject of the future investigations. 
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Branching ratio 

Predicted value 

Experiment 

Br ( D + —> 77 e + v e ) 

Br ( D + — > 77 , e + v e ) 

(14.24 ± 10.98) • 10 -4 
(1.52 ± 1.17) • 1CT 4 

(11.4 ±0.9 ±0.4) • 10 -4 [40] 
(2.16 ±0.53 ±0.07) • 10“ 4 [40] 

r (T>+ ->• 77 e + 77 e ) 

T (T> + —> r]e + v e ) 

0.10 ±0.11 

0.19 ±0.09 [41] 

Br (T)+ — > rje+i'e) 

Br (Df —> r]'e + i / e ) 

(2.40 ± 0.28)% 

(0.79 ±0.14)% 

(2.48 ± 0.29)% [41] 

(0.91 ±0.33)% [41] 

r (Z>+ -> 77 e + z 7 e ) 

T (-D+ ->• 77 e + z 7 e ) 

0.33 ±0.07 

0.36 ± 0.14[40] 

Br (B + —>• 77 e + i 7 e ) 

(0.44 ± 0.25) • 10 -4 

(0.44 ±0.23 ±0.11) • 10 -4 [42] 
(0.36 ± 0.05 ± 0.04) • 10~ 4 [43] 

Br (B + —> rj'e + v e ) 

(0.19 ±0.11) • 10 -4 

(2.66 ±0.80 ±0.56) • 10~ 4 [42] 
(0.24 ±0.08 ±0.03) • 10~ 4 [43] 

r (B + ->• 7/e + z/ e ) 

T ( B + —»■ i]e + v e ) 

0.43 ± 0.34 

0.67 ±0.24 ±0.1 [43] 

Br —> V l + l~) l=e /i 

Br (B s —> 77 r + r) 

(2.80 ± 0.36) • 10" 7 
(1.53 ± 0.18) • 10- 7 


Br(S s -^7/Z+|-) J=C)/i 

Br 

(2.85 ±0.48) • 10- 7 
(0.75 ±0.14) • 10" 7 


Br (B s —>• 77 z^P) 

Br (B s —>• 77 ' Z7p) 

(20.5 ±2.8) • 10 -7 
(14.8 ±2.0) • 10“ 7 



Table 1. Predicted branching fractions of various D^,Br s ^ —> r semileptonic decays. 
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A Explicit results for f + , f° and f J form factors at the leading order in 

B , B s —» 77W and D , D s —> 77 ^ transitions 

The leading 0(a s ) part of the /it LCSR, (3.4), has the following expression (P = 77 , 77 / 

£>(s)V K 

r = q for B q —>• P and r = s for B s —> P: for D, D s the same expressions are valid with 
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the replacement mb —> m c ): 


F 0tBr ^p(q 2 ,M 2 , sfi) = m 2 b 



2_ g 2 i l+rn? pU u ( Fp^^p(u) 
'n M - < - 


U 


1 


m b ~ <T + u 2 mj 


{vi l + q 2_ u 2 m 2 p) d&(u) 


4 um z p ml (r)a, , 

^P (“) 


rnf — q z + rrm 


+ - 


F, 


(r) 


— <r + irmf, 

2 


■#<$(«) + fi - _ 2 ^ 


m 


u 


4 mr — q 2 + rrrm 


du 2 


mg — + u z nip 

d 


6um 2 P 


+ 


dvip4p(v) 


12 amp 


m? — g 2 + u 2 m 2 p du (m? — g 2 + u 2 m 2 p ) 2 



tku) 


+ 


+ 


2um'p(ml — q 2 — u 2 m 2 P ) / d 


6umi 


(ml — q 2 + u 2 m 2 p ) 2 \du m 2 — q 2 + u 2 m 


m.fF^e f dkp(u) 




4(ra 2 — g 2 + rap) 2 \ du 


u—tl 


du 


04p(“) 


(A.l) 


where u = 1 — it, no = ^</ 2 — + mp + y/(g 2 — + mp) 2 — 4 m 2 p (q 2 — m 2 )j /(2m|>), 

= cos (j>f q /y/2,F^ = sin0/ ? /\/2, = -sin </>/,, F^ = cos 0/ s , and similarly 

for the two-particle twist-three DAs </>|p: </> 3 ^ p,cr = cos /v% = s ™tfiq /a/2, 

03 ^ P,<T = -sin^^f, 4^ P,<T = cos^^f. Also, 

= <4? w (») + . 

Jr 

tF) = <?"(») + t&“F) . ( a .2) 

Jr 


In the case of the twist-2 DA, we will express the decay constants _F 
and take the different evolution of /p 1 ^ and /p 8 * into account: 


( 9 ) 


in the SO basis 


fp\d) = fp\r o), 

1 + -¥ (“«(a 0 - «a(Mo)) 
ppo 


a)/ 


/p V) = fp Vo) 


(A.3) 


- 29 



at ?72o = 1 GeV, the energy at which the FKS parameters are determined and 

4 ,} ^2P = C q P j= + / 8 ) fo P , 


F 




4 g p = C s p-j= (/ 1 - V2f l + / 8 ) fop , 


with C q = cos0/\/2, C q , = sin^/y^, C® = — sin^>, (=, = cos</> and 


(A.4) 


02 r? = fo^' = 6 u(l - «) ( 1 + a i C i /2 ( 2u ~ !) 

*= 2,4 


(A.5) 


Numerically, 


/i(mo) = (1-17 ± 0.03)/^, 

/ 8 (/i 0 ) = (1.26 ±0.04)/*. 

The short-hand notations introduced for the integrals over three-particle DA’s are 2 * : 

U 1 

hr{u) = J dai j [4 vp -q — 3(1 - 2v)um 2 P \ 4> 3r (ai) 

0 (u-ai)/(l-ai) 

1 

i OVT/O’)/ 


«2 = 1 — a\ — « 3 , 

a 3 = (u - ai)/v 


4p ('“) = J dai J 

0 (u—ai)/(l—ai) 

1 


dp 

p 


2 *£J(ai) - ^>4p(c>!j) + 2 ^( ai ) - $S(a;) 


CK2 — 1 — CCi — a 3 , 
as = (u — a\) / v 


l[ r p(u ) = / dai 


dp 

v 


^4p( a i) + ^4p( a f) + ^4p( a i) + d^p^i) 

7 

0 (u—ai)/(l—ai) 

The leading order LCSR for f pK + f PK , (3.5), has the form 


J ! /^WP( U ) 1 d^(«) 

M \ - | - -|— 


Fo,B r ^p(q ,M , s Q ) = m b J due 

UQ 


(A-6) 

a 2 — 1 — ai — a 3 , 
a 3 = (u- ai)/v 




pOO 


+ 


m b — q 2 + u 2 m 2 p 


4tp(«) - 


2um 2 p 


m 2 ~ q 2 + u 2 m 2 p J q 
pOO \ 


2m r vib \ u 6u du 


dvip4p(v) 


old 2uTTtp \ / Fop \ ~ , * 

+mp [du m 2 -q 2 + u 2 m 2 p ) [ rnh F (P) ) 3r U 


+ - 


2um 2 p 


d 2 


6um 2 p 


d 


+ 


12u 2 rrip 


m 2 ~ q 2 + u 2 m 2 p \du 2 m 2 — q 2 + u 2 m 2 p du (m 2 — q 2 + u 2 m 2 p ) 2 


dt,lth) 


(A.7) 


2 In the paper [11], dealing with Br s \ —>• K form factors, in eq.(A.2) there was a misprint in the function 

I 3 k(u), the factor of 3 was missing. The correct expression has the same form as Isrfo) given here in (A.6). 
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where 


hr{u) = J dai J 

0 (u-ai)/(l-ai) 


dv 


— [(3 - 2v)\ $ 3 r(«i) 


a 2 = 1 — Oil — Q!3, 

«3 = (u — Ot\)/v 


Finally, the leading order LCSR for the penguin form factor, (3.6), reads: 

r m?—g 2 u+m : i > uu f F^ ( 6 ^ (u) 

= m b I due -^- P ^ 2p( } 


u 


UQ 


m b 


dcfrzp' 7 (u) 


2Uml &{u) 


6m r (m% — q 2 + u 2 m 2 p ) l du m 2 — q 2 + u 2 m 2 p 


?(r) 


+ 


2umi 


m 2 — q 2 + u 2 m 2 p \du m 2 — q 2 + u 2 m 2 p J \4 


704p(«) - A P T ( u ) 


ml u 


6 umi 


+ 


12um A P 


+ 


4(m^ — q 2 + u 2 m 2 p ) \du 2 m 2 — q 2 + u 2 m 2 p du ( m 2 — q 2 + u 2 m 2 p ) 2 


(A.8) 


04P(“) 


4(m| - g 2 + m 2 p ) 2 \ du J u ^i 


(A.9) 


and 


r ip T ( u ) = / 


du 

v 


0 (u-ai)/(l-ai) 


2*S5(oi) - (1 - 2u)$$(a i ) 


+2(1 - 2«)^5(Qi) - $2(«i) 


0^2 — 1 Ol 0:3, 

0:3 = (u — a{)/v 


(A.10) 


1 0 j 1 

The expressions for f D ’ ’ from factors follows from above, by replacing m b by m c . 


\s) 


B Parameters used in the calculation 


In this appendix we summarize the parameters used in the calculation of /m^O) form factors 
as well as in the calculation of /m decay constants, Tables 2-5. In Table 6. we summarize 
meson masses, lifetimes and vector resonances used in the calculation of phenomenological 
predictions for semileptonic M —> decays. 
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Parameter 

Value at /x = 1 GeV 

oj 

0.17 ± 0.08 [47] 

a\ 

0.06 ± 0.10 [47] 

a >4 

0 

h-K 

0.0045 ± 0.0015 GeV 2 


-1.5 ±0.7 

si 

0.18 ±0.06 GeV 2 

^47r 

0.2 ±0.1 

fzK 

0.0045 ± 0.0015 GeV 2 

^3 K 

-1.2 ±0.7 

A 3 K 

1.6 ±0.4 

Sr 

0.20 ± 0.06 GeV 2 

U)4K 

0.2 ±0.1 

H4K 

-0.09 ±0.02 


Table 2. Input parameters in DA’s [10, 21]. 


Parameter 

Value 

m b {m h ) 
m c (m c ) 
m u ( 2 GeV) 
m d (2 GeV) 
m(2 GeV) = m “+ md 
m s {2 GeV) 

4.18 ±0.03 GeV 

1.275 ± 0.025 GeV 
2.3±{j;£ MeV 
4 -8±g;| MeV 

3.5 tol MeV 

95 ± 5 MeV 

(qq)(l GeV) 
(ss)/(qq) 

{a s /ir GG) 
ml 

-(2461^ MeV) 3 
0.8 ±0.3 
0 .012+°;°°® GeV 4 
0.8 ± 0.2 GeV 2 

a s (M z ) 

0.1176 ±0.002 


Table 3. Quark masses and additional input parameters for the fB (a) and fo {s) sum rules. 
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Decay constant 

LCSR [28] 

This work 

Fitted M 2 and so (GeV 2 ) 

u 


130.7 


Ik 


155 


fn 

201±H 

191 ±9 

M 2 = 2, s 0 = 5 

fo a 

000+13 

^^o_23 

219 ± 7 

M 2 = 2, s 0 = 6.3 

fs 


215 ±7 

M 2 = 5, s 0 = 35.6 

Ib s 

242 

246 ±8 

M 2 = 5.1, s 0 = 35.5 

fp s 

Id 

1 1 r+0.04 
i - lo -0.05 

1.15 ±0.05 


JB S 

_ Ib _ 

i -| 7+0.04 
1,1 ' —0.03 

1.16 ±0.05 



Table 4. Decay constants used in the paper, the values are in MeV. The decay constants of heavy 
mesons are obtained from the two-point SR at 0(a s ) and agree with those from [10, 11, 28]. The 
quoted errors are coming only from the variation of so and the Borel parameter M, since other 
errors will cancel in the calculation of the form factors. For comparison the recent more complete 
LCSR results from [28] are given. 


Transition 
B ->• 7/W 
B s r/ 
B s —t rf 
D —> T] 

D —> rj 
D s -*r\ 
D s -f rf 


Fitted M 2 and so parameters of LCSR for / +,0,T 
M 2 b = 18 ± 2 GeV 2 , = 37 ± 0.5 GeV 2 , 
M 2 s1] = 17 ± 1 GeV 2 , Sq sV = 37.5 ± 0.5 GeV 2 
M 2 sV , = 18 ±2 GeV 2 , sJ W = 38 ± 0.5 GeV 2 
= 5.2 ± 0.8GeV 2 , s° v = 7 ± 0.2GeV 2 
M 2 Dv , = 5 ± 1 GeV 2 sJ v = 5.5 ± 0.3 GeV 2 
M 2 JsTi = 8 ± 0.2 GeV 2 , Sq sV = 7.8 ± 0.2 GeV 2 
M 2 , = 6 ± 1 GeV 2 , Sq sV ' = 7.5 ± 0.5 GeV 2 


Table 5. Fitted Borel parameters M 2 and the continuum thresholds so for each of the decays used 
to obtain the predicted form factors in the text. 


Mass 

Value (GeV) 

Resonance 

Mass value (GeV) 

Lifetime 

Value (ps) 

m B + 

5.2792 

m B *{ 1“) 

5.3252 


1.638 ±0.004 

mp Js 

5.3667 

m B *( 1“) 

5.4154 

T B S 

1.512 ±0.007 

m D + 

1.8696 

m D *(l~) 

2.0102 

1~D+ 

1.040 ±0.007 

m D+ 

1.9685 

m D *{ 1") 

2.1121 

T D+ 

0.500 ±0.007 

m n o 

0.1359 





m K o 

0.4976 





m n 

0.5478 





rrirj/ 

0.9577 






Table 6. Meson masses and lifetimes. The vector meson resonances m* H are used in the extrapo¬ 
lation formula for q 2 -dependence of the form factors (6.1). 


- 33 - 






















[4] P. Colangelo, F. De Fazio and W. Wang, “Nonleptonic B s to charmonium decays: analyses 
in pursuit of determining the weak phase /3 S ,” Phys. Rev. D 83 (2011) 094027 
[arXiv:1009.4612 [hep-ph]]. 

[5] R. Fleischer, R. Knegjens and G. Ricciardi, “Exploring CP Violation and 77 - 77 ' Mixing with 
the B ad —> J/ipri^ Systems,” Eur. Phys. J. C 71 (2011) 1798 [arXiv: 1110.5490 [hep-ph]]. 

[ 6 ] R. Aaij et al. [LHCb Collaboration], “Study of 77 — rf mixing from measurement of 
B° (s) -» J/i/h/d decay rates,” JHEP 1501 (2015) 024 [arXiv:1411.0943 [hep-ex]]. 

[7] I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, “Radiative Decay Sigma-I->■ P 

Gamma In Quantum Chromodynamics,” Nucl. Phys. B312 (1989) 509; V. M. Braun and 
I. E. Filyanov, “QCD Sum Rules In Exclusive Kinematics And Pion Wave Function,’ Z. 
Phys. C44 (1989) 157; V. L. Chernyak and I. R. Zhitnitsky, “B Meson Exclusive Decays 
Into Baryons,” Nucl. Phys. B345 (1990) 137. 

[ 8 ] P. Ball and R. Zwicky, “New results on B —»• tt, K, 77 decay formfactors from light-cone sum 
rules,” Phys. Rev. D 71 (2005) 014015 [hep-pli/0406232]. 

[9] A. Khodjamirian, C. Klein, T. Mannel and N. Offen, “Semileptonic charm decays 

D —> nlu(l ) and D —> Klv{l ) from QCD Light-Cone Sum Rules,” Phys. Rev. D 80 (2009) 
114005 [arXiv:0907.2842 [hep-ph]]. 

[10] G. Duplancic, A. Khodjamirian, T. Mannel, B. Melic and N. Offen, “Light-cone sum rules 
for B —» tt form factors revisited,” JHEP 0804 (2008) 014. [arXiv:0801.1796 [hep-ph]]. 

[11] G. Duplancic and B. Melic, “J5, B(s) K form factors: An Update of light-cone sum rule 
results,” Phys. Rev. D 78 (2008) 054015 [arXiv:0805.4170 [hep-ph]]. 

[12] S. S. Agaev, V. M. Braun, N. Offen, F. A. Porkert and A. Schfer, “Transition form factors 

7*7 ->• 77 and 7*7 77 ' in QCD,” Phys. Rev. D 90 (2014) 7, 074019 [arXiv:1409.4311 

[hep-ph]]. 

[13] P. Ball and G. W. Jones, “B -» 77 W Form Factors in QCD,” JHEP 0708 (2007) 025 
[arXiv:0706.3628 [hep-ph]]. 

[14] P. Kroll and K. Passek-Kumericki, “The Two gluon components of the eta and eta-prime 
mesons to leading twist accuracy,” Phys. Rev. D 67 (2003) 054017 [hep-ph/0210045]. 

[15] T. Feldmann, P. Kroll and B. Stech, “Mixing and decay constants of pseudoscalar mesons,” 
Phys. Rev. D 58 (1998) 114006 [hep-ph/9802409]. 

[16] F. Ambrosino, A. Antonelli, M. Antonelli, F. Archilli, P. Beltrame, G. Bencivenni, 

S. Bertolucci and C. Bini et al., “A Global fit to determine the pseudoscalar mixing angle 
and the gluonium content of the eta-prime meson,” JHEP 0907 (2009) 105 
[arXiv:0906.3819 [hep-ph]]. 

[17] C. Di Donato, G. Ricciardi and I. Bigi, “77 — 77 ' Mixing - From electromagnetic transitions to 
weak decays of charm and beauty hadrons,” Phys. Rev. D 85 (2012) 013016 
[arXiv:1105.3557 [hep-ph]]. 

[18] C. Michael et al. [ETM Collaboration], “77 and 77 ' mixing from Lattice QCD,” Phys. Rev. 
Lett. Ill (2013) 18, 181602 [arXiv:1310.1207 [hep-lat]]. 

[19] V. N. Baier and A. G. Grozin, “Meson Wave Functions With Two Gluon States,” Nucl. 
Phys. B 192 (1981) 476. 


- 34 - 



[20] B. Melic, B. Nizic and K. Passek, “A Note on the factorization scale independence of the 
PQCD predictions for exclusive processes,” Eur. Phys. J. C 36 (2004) 453 [hep-ph/0107311]. 

[21] P. Ball, V. M. Braun and A. Lenz, “Higher-twist distribution amplitudes of the K meson in 
QCD,” JHEP 0605 (2006) 004. 

[22] M. Beneke and M. Neubert, “Flavor singlet B decay amplitudes in QCD factorization,” 
Nucl. Phys. B 651 (2003) 225 [hep-ph/0210085], 

[23] M. Jamin and B. O. Lange, “f(B) and f(B/s) from QCD sum rules,” Phys. Rev. D 65, 
056005 (2002). 

[24] G. P. Lepage and S. J. Brodsky, “Exclusive Processes In Quantum Chromodynamics: 
Evolution Equations For Hadronic Wave Functions And The Form-Factors Of Mesons,” 
Phys. Lett. B 87, 359 (1979); “Exclusive Processes In Perturbative Quantum 
Chromodynamics,” Phys. Rev. D 22, 2157 (1980); A. V. Efremov and A. V. Radyushkin, 
“Factorization And Asymptotical Behavior Of Pion Form-Factor In QCD,” Phys. Lett. B 
94, 245 (1980); “Asymptotical Behavior Of Pion Electromagnetic Form-Factor In QCD,” 
Theor. Math. Phys. 42, 97 (1980). 

[25] M. V. Terentev, “Factorization In Exclusive Processes. Form-factor Of Singlet Mesons In 
Quantum Chromodynamics,” Sov. J. Nucl. Phys. 33 (1981) 911 [Yad. Fiz. 33 (1981) 1692]. 

[26] B. Melic, B. Nizic and K. Passek, “BLM scale setting for the pion transition form-factor,” 
Phys. Rev. D 65 (2002) 053020 [hep-ph/0107295]. 

[27] A. Khodjamirian, R. Riickl, S. Weinzierl and O. I. Yakovlev, “Perturbative QCD correction 
to the B —» 7r transition form factor,” Phys. Lett. B 410 (1997) 275. 

[28] P. Gelhausen, A. Khodjamirian, A. A. Pivovarov and D. Rosenthal, “Decay constants of 
heavy-light vector mesons from QCD sum rules,” Phys. Rev. D 88 (2013) 1, 014015 
[Erratum-ibid. D 89 (2014) 9, 099901] [arXiv: 1305.5432 [hep-ph]]. 

[29] R. J. Dowdall et al. [HPQCD Collaboration], “B-Meson Decay Constants from Improved 
Lattice Nonrelativistic QCD with Physical u, d, s, and c Quarks,” Phys. Rev. Lett. 110 
(2013) 22, 222003 [arXiv:1302.2644 [hep-lat]]. A. Bazavov et al. [Fermilab Lattice and 
MILC Collaborations], “Charmed and light pseudoscalar meson decay constants from 
four-flavor lattice QCD with physical light quarks,” Phys. Rev. D 90 (2014) 7, 074509 
[arXiv:1407.3772 [hep-lat]]. N. Carrasco, P. Dimopoulos, R. Frezzotti, P. Larni, V. Lubicz, 

F. Nazzaro, E. Picca and L. Riggio et al., “Leptonic decay constants fK, fD and fDs with Nf 
= 2+1+1 twisted-mass lattice QCD,” arXiv: 1411.7908 [hep-lat]. 

[30] N. Offen, F. A. Porkert and A. Schfer, “Light-cone sum rules for the D( s ) -+ rj^’lui form 
factor,” Phys. Rev. D 88 (2013) 3, 034023 [arXiv: 1307.2797], 

[31] Y. Y. Charng, T. Kurimoto and H. n. Li, “Gluonic contribution to B — > rj^ form factors,” 
Phys. Rev. D 74 (2006) 074024 [Erratum-ibid. D 78 (2008) 059901] [hep-ph/0609165]. 

[32] T. M. Aliev, I. Kanik and A. Ozpineci, “Semileptonic B -+ i)lv decay in light cone QCD,” 
Phys. Rev. D 67 (2003) 094009 [hep-ph/0210403]. 

[33] T. M. Aliev and M. Savci, “Form-factor for penguin-induced B —» r] transitions in light cone 
QCD sum rules,” Eur. Phys. J. C 29 (2003) 515 [hep-ph/0302187], 

[34] P. Colangelo and F. De Fazio, “D(s) decays to eta and eta-prime final states: A 
Phenomenological analysis,” Phys. Lett. B 520 (2001) 78 [hep-ph/0107137]. 


- 35 - 



[35] K. Azizi, R. Khosravi and F. Falaliati, “Exclusive D s —> ( 77 , rf)lv decays in light cone 
QCD,” J. Phys. G 38 (2011) 095001 [arXiv: 1011.6046 [hep-ph]]. 

[36] G. S. Bali, S. Collins, S. Durr and I. Kanamori, U D S —► r), rf semileptonic decay form factors 
with disconnected quark loop contributions,” Phys. Rev. D 91 (2015) 014503 
[arXiv:1406.5449 [hep-lat]]. 

[37] D. Becirevic and A. B. Kaidalov, “Comment on the heavy —> light form-factors,” Phys. 
Lett. B 478 (2000) 417 [hep-ph/9904490], 

[38] G. Buchalla, A. J. Buras and M. E. Lautenbacher, “Weak decays beyond leading 
logarithms,” Rev. Mod. Phys. 68 (1996) 1125 [hep-ph/9512380]. 

[39] M. V. Carlucci, P. Colangelo and F. De Fazio, “Rare B(s) decays to eta and eta-prime final 
states,” Phys. Rev. D 80 (2009) 055023 [arXiv:0907.2160 [hep-ph]]. 

[40] J. Yelton et al. [CLEO Collaboration], “Studies of D + —y (if, 77 , cf))e + v e ” , Phys. Rev. D 84 
(2011) 032001 [arXiv:1011.1195 [hep-ex]]. 

[41] J. Yelton et al. [CLEO Collaboration], “Absolute Branching Fraction Measurements for 
Exclusive D(s) Semileptonic Decays,” Phys. Rev. D 80 (2009) 052007 [arXiv:0903.0601 
[hep-ex]]. 

[42] N. E. Adam et al. [CLEO Collaboration], “A Study of Exclusive Charmless Semileptonic B 
Decay and \V(ub)\”, Phys. Rev. Lett. 99 (2007) 041802 [hep-ex/0703041 [HEP-EX]]. 

[43] P. del Amo Sanchez et al. [BaBar Collaboration], “Measurement of the B° —» n e £ + v and 
B + —> rf ^>l + v Branching Fractions, the B° —> n~i + v and B + —> t]£ + i/ Form-Factor Shapes, 
and Determination of |14b|,” Phys. Rev. D 83 (2011) 052011 [arXiv:1010.0987 [hep-ex]]. 

[44] K. Azizi, R. Khosravi and F. Falahati, “Rare Semileptonic B s Decays to 77 and if mesons in 
QCD,” Phys. Rev. D 82 (2010) 116001 [arXiv: 1008.3175 [hep-ph]]. 

[45] H. M. Choi, “Exclusive Rare B s —» (K, 77, r)')l + £~ Decays in the Light-Front Quark Model,” 
J. Phys. G 37 (2010) 085005 [arXiv:1002.0721 [hep-ph]]. 

[46] R. N. Faustov and V. O. Galkin, “Rare B s decays in the relativistic quark model,” Eur. 
Phys. J. C 73 (2013) 10, 2593 [arXiv:1309.2160 [hep-ph]]. 

[47] A. Khodjamirian, T. Mannel, N. Offen and Y.-M. Wang, a B 1 t&i Width and \ V u f\ from 

QCD Light-Cone Sum Rules,” Phys. Rev. D 83 (2011) 094031 [arXiv:1103.2655 [hep-ph]]. 


- 36 - 



